(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[    118938,       3653]
NotebookOptionsPosition[    116365,       3568]
NotebookOutlinePosition[    116703,       3583]
CellTagsIndexPosition[    116660,       3580]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{"\[Rho]", " ", "=", " ", 
  RowBox[{"U", "[", "1", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"u", " ", "=", " ", 
  RowBox[{
   RowBox[{"U", "[", "2", "]"}], "/", 
   RowBox[{"U", "[", "1", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"v", " ", "=", " ", 
  RowBox[{
   RowBox[{"U", "[", "3", "]"}], "/", 
   RowBox[{"U", "[", "1", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"e", " ", "=", " ", 
  RowBox[{"U", "[", "4", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"p", " ", "=", " ", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"\[Gamma]", "-", "1"}], ")"}], 
   RowBox[{"(", 
    RowBox[{"e", "-", 
     RowBox[{"\[Rho]", " ", 
      RowBox[{
       RowBox[{"u", "^", "2"}], "/", "2"}]}], "-", 
     RowBox[{"\[Rho]", " ", 
      RowBox[{
       RowBox[{"v", "^", "2"}], "/", "2"}]}]}], 
    ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"F", "=", 
  RowBox[{"{", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", " ", "u"}], ",", 
     RowBox[{
      RowBox[{"\[Rho]", " ", 
       RowBox[{
        RowBox[{"(", "u", ")"}], "^", "2"}]}], "+", "p"}], ",", 
     RowBox[{"\[Rho]", " ", "u", " ", "v"}], ",", 
     RowBox[{"u", 
      RowBox[{"(", 
       RowBox[{"e", "+", "p"}], ")"}]}]}], "}"}], 
   "}"}]}], "\[IndentingNewLine]", 
 RowBox[{"G", " ", "=", 
  RowBox[{"{", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", " ", "v"}], ",", 
     RowBox[{"\[Rho]", " ", "u", " ", "v"}], ",", 
     RowBox[{
      RowBox[{"\[Rho]", " ", 
       RowBox[{"v", "^", "2"}]}], "+", "p"}], ",", 
     RowBox[{"v", 
      RowBox[{"(", 
       RowBox[{"e", "+", "p"}], ")"}]}]}], "}"}], "}"}]}]}], "Input",
 CellChangeTimes->{{3.513759174050398*^9, 3.513759474881263*^9}, {
  3.513759568378648*^9, 3.513759582101974*^9}, {3.513760158085368*^9, 
  3.513760167042832*^9}, {3.513769981211738*^9, 3.513769983468673*^9}, {
  3.513820004982169*^9, 3.513820028711392*^9}}],

Cell[BoxData[
 RowBox[{"U", "[", "1", "]"}]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403226286*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{"U", "[", "2", "]"}], 
  RowBox[{"U", "[", "1", "]"}]]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403230575*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{"U", "[", "3", "]"}], 
  RowBox[{"U", "[", "1", "]"}]]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403233953*^9}],

Cell[BoxData[
 RowBox[{"U", "[", "4", "]"}]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403237073*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", 
     FractionBox[
      SuperscriptBox[
       RowBox[{"U", "[", "2", "]"}], "2"], 
      RowBox[{"2", " ", 
       RowBox[{"U", "[", "1", "]"}]}]]}], "-", 
    FractionBox[
     SuperscriptBox[
      RowBox[{"U", "[", "3", "]"}], "2"], 
     RowBox[{"2", " ", 
      RowBox[{"U", "[", "1", "]"}]}]], "+", 
    RowBox[{"U", "[", "4", "]"}]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403240978*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"U", "[", "2", "]"}], ",", 
    RowBox[{
     FractionBox[
      SuperscriptBox[
       RowBox[{"U", "[", "2", "]"}], "2"], 
      RowBox[{"U", "[", "1", "]"}]], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         FractionBox[
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"], 
          RowBox[{"2", " ", 
           RowBox[{"U", "[", "1", "]"}]}]]}], "-", 
        FractionBox[
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"], 
         RowBox[{"2", " ", 
          RowBox[{"U", "[", "1", "]"}]}]], "+", 
        RowBox[{"U", "[", "4", "]"}]}], ")"}]}]}], ",", 
    FractionBox[
     RowBox[{
      RowBox[{"U", "[", "2", "]"}], " ", 
      RowBox[{"U", "[", "3", "]"}]}], 
     RowBox[{"U", "[", "1", "]"}]], ",", 
    FractionBox[
     RowBox[{
      RowBox[{"U", "[", "2", "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"U", "[", "4", "]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"U", "[", "2", "]"}], "2"], 
             RowBox[{"2", " ", 
              RowBox[{"U", "[", "1", "]"}]}]]}], "-", 
           FractionBox[
            SuperscriptBox[
             RowBox[{"U", "[", "3", "]"}], "2"], 
            RowBox[{"2", " ", 
             RowBox[{"U", "[", "1", "]"}]}]], "+", 
           RowBox[{"U", "[", "4", "]"}]}], ")"}]}]}], ")"}]}], 
     RowBox[{"U", "[", "1", "]"}]]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403246357*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"U", "[", "3", "]"}], ",", 
    FractionBox[
     RowBox[{
      RowBox[{"U", "[", "2", "]"}], " ", 
      RowBox[{"U", "[", "3", "]"}]}], 
     RowBox[{"U", "[", "1", "]"}]], ",", 
    RowBox[{
     FractionBox[
      SuperscriptBox[
       RowBox[{"U", "[", "3", "]"}], "2"], 
      RowBox[{"U", "[", "1", "]"}]], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         FractionBox[
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"], 
          RowBox[{"2", " ", 
           RowBox[{"U", "[", "1", "]"}]}]]}], "-", 
        FractionBox[
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"], 
         RowBox[{"2", " ", 
          RowBox[{"U", "[", "1", "]"}]}]], "+", 
        RowBox[{"U", "[", "4", "]"}]}], ")"}]}]}], ",", 
    FractionBox[
     RowBox[{
      RowBox[{"U", "[", "3", "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"U", "[", "4", "]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"U", "[", "2", "]"}], "2"], 
             RowBox[{"2", " ", 
              RowBox[{"U", "[", "1", "]"}]}]]}], "-", 
           FractionBox[
            SuperscriptBox[
             RowBox[{"U", "[", "3", "]"}], "2"], 
            RowBox[{"2", " ", 
             RowBox[{"U", "[", "1", "]"}]}]], "+", 
           RowBox[{"U", "[", "4", "]"}]}], ")"}]}]}], ")"}]}], 
     RowBox[{"U", "[", "1", "]"}]]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.513759179265598*^9, 3.513759207300768*^9}, {
   3.513759237784863*^9, 3.513759244272892*^9}, 3.513759301389891*^9, {
   3.513759393293308*^9, 3.513759422036931*^9}, 3.513759477135175*^9, 
   3.513759582674105*^9, 3.513760167713748*^9, 3.513767219030814*^9, 
   3.513767447093402*^9, 3.513769941477512*^9, 3.513813403251489*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"A", "=", 
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"F", ",", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"U", "[", "1", "]"}], ",", 
        RowBox[{"U", "[", "2", "]"}], ",", 
        RowBox[{"U", "[", "3", "]"}], ",", 
        RowBox[{"U", "[", "4", "]"}]}], "}"}], "}"}]}], "]"}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"B", "=", 
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"G", ",", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"U", "[", "1", "]"}], ",", 
        RowBox[{"U", "[", "2", "]"}], ",", 
        RowBox[{"U", "[", "3", "]"}], ",", 
        RowBox[{"U", "[", "4", "]"}]}], "}"}], "}"}]}], "]"}], "//", 
   "Simplify"}]}]}], "Input",
 CellChangeTimes->{{3.513759479893178*^9, 3.513759504142008*^9}, {
  3.513759536865459*^9, 3.513759550713114*^9}, {3.513759617146509*^9, 
  3.513759617428744*^9}, {3.513759741530364*^9, 3.51375974745009*^9}, {
  3.513760436892497*^9, 3.51376049543132*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"U", "[", "3", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", 
      FractionBox[
       RowBox[{"U", "[", "2", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
          RowBox[{"\[Gamma]", " ", 
           RowBox[{"U", "[", "1", "]"}], " ", 
           RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"3", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"2", " ", "\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"\[Gamma]", " ", 
        RowBox[{"U", "[", "2", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.513813529456325*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", "0", ",", "1", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"U", "[", "3", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", 
      FractionBox[
       RowBox[{"U", "[", "2", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "3", "]"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
          RowBox[{"\[Gamma]", " ", 
           RowBox[{"U", "[", "1", "]"}], " ", 
           RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{"3", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"2", " ", "\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
      FractionBox[
       RowBox[{"\[Gamma]", " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.513813529752999*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"A", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"U", "[", "3", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", 
      FractionBox[
       RowBox[{"U", "[", "2", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
          RowBox[{"\[Gamma]", " ", 
           RowBox[{"U", "[", "1", "]"}], " ", 
           RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"3", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"2", " ", "\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"\[Gamma]", " ", 
        RowBox[{"U", "[", "2", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.513770176145311*^9, 3.513770177230981*^9}, {
  3.51377023897528*^9, 3.513770248088954*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
      RowBox[{"2", " ", 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "2", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], ",", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
     FractionBox[
      RowBox[{"U", "[", "3", "]"}], 
      RowBox[{"U", "[", "1", "]"}]], ",", 
     FractionBox[
      RowBox[{"U", "[", "2", "]"}], 
      RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"U", "[", "2", "]"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
      SuperscriptBox[
       RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"3", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
        RowBox[{"2", " ", "\[Gamma]", " ", 
         RowBox[{"U", "[", "1", "]"}], " ", 
         RowBox[{"U", "[", "4", "]"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
     FractionBox[
      RowBox[{"\[Gamma]", " ", 
       RowBox[{"U", "[", "2", "]"}]}], 
      RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513770204971247*^9, 3.513770255443044*^9, 
  3.513813585720373*^9, 3.513814832357432*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", "0", ",", "1", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      FractionBox[
       RowBox[{"U", "[", "3", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", 
      FractionBox[
       RowBox[{"U", "[", "2", "]"}], 
       RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        RowBox[{"U", "[", "1", "]"}]]}], ",", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "3", "]"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           SuperscriptBox[
            RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
          RowBox[{"\[Gamma]", " ", 
           RowBox[{"U", "[", "1", "]"}], " ", 
           RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         RowBox[{"U", "[", "2", "]"}], " ", 
         RowBox[{"U", "[", "3", "]"}]}], 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{"3", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"2", " ", "\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
      FractionBox[
       RowBox[{"\[Gamma]", " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.513770179672063*^9, 3.51377018019501*^9}, {
  3.513770250682831*^9, 3.513770253131069*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "1", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
     FractionBox[
      RowBox[{"U", "[", "3", "]"}], 
      RowBox[{"U", "[", "1", "]"}]], ",", 
     FractionBox[
      RowBox[{"U", "[", "2", "]"}], 
      RowBox[{"U", "[", "1", "]"}]], ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "3", "]"}], "2"]}]}], 
      RowBox[{"2", " ", 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]}]], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "2", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       RowBox[{"U", "[", "1", "]"}]]}], ",", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"U", "[", "3", "]"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          SuperscriptBox[
           RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
         RowBox[{"\[Gamma]", " ", 
          RowBox[{"U", "[", "1", "]"}], " ", 
          RowBox[{"U", "[", "4", "]"}]}]}], ")"}]}], 
      SuperscriptBox[
       RowBox[{"U", "[", "1", "]"}], "3"]], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
        RowBox[{"U", "[", "2", "]"}], " ", 
        RowBox[{"U", "[", "3", "]"}]}], 
       SuperscriptBox[
        RowBox[{"U", "[", "1", "]"}], "2"]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "2", "]"}], "2"]}], "+", 
        RowBox[{"3", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
         SuperscriptBox[
          RowBox[{"U", "[", "3", "]"}], "2"]}], "-", 
        RowBox[{"2", " ", "\[Gamma]", " ", 
         RowBox[{"U", "[", "1", "]"}], " ", 
         RowBox[{"U", "[", "4", "]"}]}]}], 
       RowBox[{"2", " ", 
        SuperscriptBox[
         RowBox[{"U", "[", "1", "]"}], "2"]}]]}], ",", 
     FractionBox[
      RowBox[{"\[Gamma]", " ", 
       RowBox[{"U", "[", "3", "]"}]}], 
      RowBox[{"U", "[", "1", "]"}]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513770202173314*^9, 3.513770257879234*^9, 
  3.51381358962136*^9, 3.513814835541728*^9}]
}, Open  ]],

Cell[BoxData["Quit"], "Input",
 CellChangeTimes->{{3.513760746375796*^9, 3.513760746816289*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"A", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"U", "[", "1", "]"}], "\[Rule]", "\[Rho]"}], ",", 
      RowBox[{
       RowBox[{"U", "[", "2", "]"}], "\[Rule]", 
       RowBox[{"\[Rho]", " ", "u"}]}], ",", " ", 
      RowBox[{
       RowBox[{"U", "[", "3", "]"}], "\[Rule]", " ", 
       RowBox[{"\[Rho]", " ", "v"}]}], ",", " ", 
      RowBox[{
       RowBox[{"U", "[", "4", "]"}], "\[Rule]", "e"}]}], "}"}]}], "//", 
   "Simplify"}], "//", "MatrixForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"B", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"U", "[", "1", "]"}], "\[Rule]", "\[Rho]"}], ",", 
      RowBox[{
       RowBox[{"U", "[", "2", "]"}], "\[Rule]", 
       RowBox[{"\[Rho]", " ", "u"}]}], ",", " ", 
      RowBox[{
       RowBox[{"U", "[", "3", "]"}], "\[Rule]", " ", 
       RowBox[{"\[Rho]", " ", "v"}]}], ",", " ", 
      RowBox[{
       RowBox[{"U", "[", "4", "]"}], "\[Rule]", "e"}]}], "}"}]}], "//", 
   "Simplify"}], "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.513760565070102*^9, 3.513760631712487*^9}, {
  3.513760770634752*^9, 3.513760790361951*^9}, {3.51376093745351*^9, 
  3.513760946892463*^9}, {3.513770209407828*^9, 3.513770210638293*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "1", "0", "0"},
     {
      RowBox[{
       FractionBox["1", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SuperscriptBox["u", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
         RowBox[{
          SuperscriptBox["v", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], 
      RowBox[{
       RowBox[{"-", "u"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], 
      RowBox[{"v", "-", 
       RowBox[{"v", " ", "\[Gamma]"}]}], 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]},
     {
      RowBox[{
       RowBox[{"-", "u"}], " ", "v"}], "v", "u", "0"},
     {
      FractionBox[
       RowBox[{"u", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             SuperscriptBox["u", "2"], "+", 
             SuperscriptBox["v", "2"]}], ")"}], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
         ")"}]}], "\[Rho]"], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"3", " ", 
            SuperscriptBox["u", "2"]}], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
       RowBox[{"2", " ", "\[Rho]"}]], 
      RowBox[{
       RowBox[{"-", "u"}], " ", "v", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
      RowBox[{"u", " ", "\[Gamma]"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.513813597084156*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "1", "0"},
     {
      RowBox[{
       RowBox[{"-", "u"}], " ", "v"}], "v", "u", "0"},
     {
      RowBox[{
       FractionBox["1", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SuperscriptBox["v", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
         RowBox[{
          SuperscriptBox["u", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], 
      RowBox[{"u", "-", 
       RowBox[{"u", " ", "\[Gamma]"}]}], 
      RowBox[{
       RowBox[{"-", "v"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}]},
     {
      FractionBox[
       RowBox[{"v", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             SuperscriptBox["u", "2"], "+", 
             SuperscriptBox["v", "2"]}], ")"}], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
         ")"}]}], "\[Rho]"], 
      RowBox[{
       RowBox[{"-", "u"}], " ", "v", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           RowBox[{"3", " ", 
            SuperscriptBox["v", "2"]}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
       RowBox[{"2", " ", "\[Rho]"}]], 
      RowBox[{"v", " ", "\[Gamma]"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.513813597302476*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 TagBox[
  RowBox[{"A", "=", 
   RowBox[{"(", "\[NoBreak]", GridBox[{
      {"0", "1", "0", "0"},
      {
       RowBox[{
        FractionBox["1", "2"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           SuperscriptBox["u", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
          RowBox[{
           SuperscriptBox["v", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], 
       RowBox[{
        RowBox[{"-", "u"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], 
       RowBox[{"v", "-", 
        RowBox[{"v", " ", "\[Gamma]"}]}], 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}]},
      {
       RowBox[{
        RowBox[{"-", "u"}], " ", "v"}], "v", "u", "0"},
      {
       FractionBox[
        RowBox[{"u", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}],
           ")"}]}], "\[Rho]"], 
       FractionBox[
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"3", " ", 
             SuperscriptBox["u", "2"]}], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
        RowBox[{"2", " ", "\[Rho]"}]], 
       RowBox[{
        RowBox[{"-", "u"}], " ", "v", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       RowBox[{"u", " ", "\[Gamma]"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.7]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Input",
 CellChangeTimes->{{3.513770398176755*^9, 3.513770398956717*^9}, {
   3.513813754657668*^9, 3.513813799317098*^9}, {3.51381388927334*^9, 
   3.513813907039545*^9}, 3.513813972052398*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SuperscriptBox["u", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
        RowBox[{
         SuperscriptBox["v", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], ",", 
     RowBox[{
      RowBox[{"-", "u"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], ",", 
     RowBox[{"v", "-", 
      RowBox[{"v", " ", "\[Gamma]"}]}], ",", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "u"}], " ", "v"}], ",", "v", ",", "u", ",", "0"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{"u", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
        ")"}]}], "\[Rho]"], ",", 
     FractionBox[
      RowBox[{
       RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"3", " ", 
           SuperscriptBox["u", "2"]}], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
      RowBox[{"2", " ", "\[Rho]"}]], ",", 
     RowBox[{
      RowBox[{"-", "u"}], " ", "v", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], ",", 
     RowBox[{"u", " ", "\[Gamma]"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513813977861449*^9, 3.513814845141852*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"/.", 
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"p", " ", "==", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Gamma]", "-", "1"}], ")"}], 
       RowBox[{"(", 
        RowBox[{"e", "-", 
         RowBox[{"\[Rho]", " ", 
          RowBox[{
           RowBox[{"u", "^", "2"}], "/", "2"}]}], "-", 
         RowBox[{"\[Rho]", " ", 
          RowBox[{
           RowBox[{"v", "^", "2"}], "/", "2"}]}]}], ")"}]}]}], ",", "e"}], 
    "]"}]}], "//", "FullSimplify"}]], "Input",
 CellChangeTimes->{3.513813975623715*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B", "=", 
  TagBox[
   RowBox[{"(", "\[NoBreak]", GridBox[{
      {"0", "0", "1", "0"},
      {
       RowBox[{
        RowBox[{"-", "u"}], " ", "v"}], "v", "u", "0"},
      {
       RowBox[{
        FractionBox["1", "2"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           SuperscriptBox["v", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
          RowBox[{
           SuperscriptBox["u", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], 
       RowBox[{"u", "-", 
        RowBox[{"u", " ", "\[Gamma]"}]}], 
       RowBox[{
        RowBox[{"-", "v"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}]},
      {
       FractionBox[
        RowBox[{"v", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}],
           ")"}]}], "\[Rho]"], 
       RowBox[{
        RowBox[{"-", "u"}], " ", "v", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], 
       FractionBox[
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            RowBox[{"3", " ", 
             SuperscriptBox["v", "2"]}]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
        RowBox[{"2", " ", "\[Rho]"}]], 
       RowBox[{"v", " ", "\[Gamma]"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.7]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
   Function[BoxForm`e$, 
    MatrixForm[BoxForm`e$]]]}]], "Input",
 CellChangeTimes->{{3.51377040166882*^9, 3.513770404616779*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "1", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "u"}], " ", "v"}], ",", "v", ",", "u", ",", "0"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         SuperscriptBox["v", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], "+", 
        RowBox[{
         SuperscriptBox["u", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}]}], ")"}]}], ",", 
     RowBox[{"u", "-", 
      RowBox[{"u", " ", "\[Gamma]"}]}], ",", 
     RowBox[{
      RowBox[{"-", "v"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}]}], ",", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{"v", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "e"}], " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
        ")"}]}], "\[Rho]"], ",", 
     RowBox[{
      RowBox[{"-", "u"}], " ", "v", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], ",", 
     FractionBox[
      RowBox[{
       RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          RowBox[{"3", " ", 
           SuperscriptBox["v", "2"]}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}]}], 
      RowBox[{"2", " ", "\[Rho]"}]], ",", 
     RowBox[{"v", " ", "\[Gamma]"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513770408872872*^9, 3.513813624916744*^9, 
  3.513813691421933*^9, 3.513814849142435*^9, 3.513819691643518*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"p", " ", "=", " ", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"\[Gamma]", "-", "1"}], ")"}], 
   RowBox[{"(", 
    RowBox[{"e", "-", 
     RowBox[{"\[Rho]", " ", 
      RowBox[{
       RowBox[{"u", "^", "2"}], "/", "2"}]}], "-", 
     RowBox[{"\[Rho]", " ", 
      RowBox[{
       RowBox[{"v", "^", "2"}], "/", "2"}]}]}], 
    ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"H", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"e", "+", "p"}], ")"}], "/", 
   "\[Rho]"}]}], "\[IndentingNewLine]", 
 RowBox[{"a", "=", " ", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"\[Gamma]", "  ", 
    RowBox[{"p", "/", "\[Rho]"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.513760688752114*^9, 3.513760704169449*^9}, {
  3.513761491094468*^9, 3.513761493368145*^9}, {3.513761532497335*^9, 
  3.51376153377159*^9}, {3.513767628370125*^9, 3.513767636241204*^9}, {
  3.513767681107127*^9, 3.513767685202242*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{"e", "-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]], "Output",
 CellChangeTimes->{3.513820046291673*^9, 3.513820563456255*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{"e", "+", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{"e", "-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}]}], 
  "\[Rho]"]], "Output",
 CellChangeTimes->{3.513820046291673*^9, 3.513820563463319*^9}],

Cell[BoxData[
 SqrtBox[
  FractionBox[
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
    RowBox[{"(", 
     RowBox[{"e", "-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
   "\[Rho]"]]], "Output",
 CellChangeTimes->{3.513820046291673*^9, 3.513820563469756*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Eigenvectors", "[", "B", "]"}], "//", "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.513819744929395*^9, 3.513819780914889*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"ebt", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", 
      FractionBox["1", "u"], ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      FractionBox["1", "v"], ",", 
      RowBox[{
       FractionBox["u", 
        RowBox[{"2", " ", "v"}]], "-", 
       FractionBox["v", 
        RowBox[{"2", " ", "u"}]]}], ",", "1", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{"2", " ", "\[Rho]"}], 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "+", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"2", " ", "u", " ", "\[Rho]"}], 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "+", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}]]}], ",", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}],
            "+", 
           RowBox[{
            SuperscriptBox["u", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "v", " ", 
               RowBox[{"(", 
                RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
              RowBox[{
               SqrtBox["2"], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"-", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                 "\[Gamma]", " ", "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}], "+", 
           RowBox[{
            SuperscriptBox["v", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "v", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "1"}], "+", 
                 SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], 
              "-", 
              RowBox[{
               SqrtBox["2"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"-", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                 "\[Gamma]", " ", "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"u", "-", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{"u", "+", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "+", 
            RowBox[{
             SqrtBox["2"], " ", "v", " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"-", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
               "\[Gamma]", " ", "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
          
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["u", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "+", 
            RowBox[{"v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "v"}], " ", 
                RowBox[{"(", 
                 RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
               
               RowBox[{"2", " ", 
                SqrtBox["2"], " ", 
                SqrtBox[
                 RowBox[{
                  RowBox[{"-", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                  "\[Gamma]", " ", "\[Rho]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                   ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ",", "1"}], 
     "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{"2", " ", "\[Rho]"}], 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "-", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}]]}], ",", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"2", " ", "u", " ", "\[Rho]"}], 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "-", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}]]}], ",", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}],
            "+", 
           RowBox[{
            SuperscriptBox["u", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "v", " ", 
               RowBox[{"(", 
                RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
              RowBox[{
               SqrtBox["2"], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"-", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                 "\[Gamma]", " ", "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}], "+", 
           RowBox[{
            SuperscriptBox["v", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "v", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "1"}], "+", 
                 SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], 
              "+", 
              RowBox[{
               SqrtBox["2"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"-", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                 "\[Gamma]", " ", "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"u", "-", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{"u", "+", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "-", 
            RowBox[{
             SqrtBox["2"], " ", "v", " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"-", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
               "\[Gamma]", " ", "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
          
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["u", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "-", 
            RowBox[{"v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"v", " ", 
                RowBox[{"(", 
                 RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
               
               RowBox[{"2", " ", 
                SqrtBox["2"], " ", 
                SqrtBox[
                 RowBox[{
                  RowBox[{"-", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
                  "\[Gamma]", " ", "\[Rho]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                   ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ",", "1"}], 
     "}"}]}], "}"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ebt", "[", 
   RowBox[{"[", "1", "]"}], "]"}], "=", 
  RowBox[{
   RowBox[{"ebt", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "/", 
   RowBox[{"ebt", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ebt", "[", 
   RowBox[{"[", "2", "]"}], "]"}], "=", 
  RowBox[{
   RowBox[{"ebt", "[", 
    RowBox[{"[", "2", "]"}], "]"}], "/", 
   RowBox[{"ebt", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ebt", "[", 
   RowBox[{"[", "3", "]"}], "]"}], "=", 
  RowBox[{
   RowBox[{"ebt", "[", 
    RowBox[{"[", "3", "]"}], "]"}], "/", 
   RowBox[{"ebt", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ebt", "[", 
    RowBox[{"[", "4", "]"}], "]"}], "=", 
   RowBox[{
    RowBox[{"ebt", "[", 
     RowBox[{"[", "4", "]"}], "]"}], "/", 
    RowBox[{"ebt", "[", 
     RowBox[{"[", 
      RowBox[{"4", ",", "1"}], "]"}], "]"}]}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.513819809121865*^9, 3.513819810810877*^9}, {
  3.513819886022199*^9, 3.513819890376453*^9}, {3.513820246998568*^9, 
  3.513820359249084*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", 
     FractionBox["1", "u"], ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "v"], ",", 
     RowBox[{
      FractionBox["u", 
       RowBox[{"2", " ", "v"}]], "-", 
      FractionBox["v", 
       RowBox[{"2", " ", "u"}]]}], ",", "1", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{"2", " ", "\[Rho]"}], 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
           "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{"2", " ", "u", " ", "\[Rho]"}], 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
           "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}]]}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Rho]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}], 
          "+", 
          RowBox[{
           SuperscriptBox["u", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "v", " ", 
              RowBox[{"(", 
               RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
             RowBox[{
              SqrtBox["2"], " ", 
              SqrtBox[
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
                "\[Rho]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "e"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                 ")"}]}]]}]}], ")"}]}], "+", 
          RowBox[{
           SuperscriptBox["v", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "v", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "1"}], "+", 
                SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], 
             "-", 
             RowBox[{
              SqrtBox["2"], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
              SqrtBox[
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
                "\[Rho]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "e"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                 ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"u", "-", "v"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{"u", "+", "v"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
           "+", 
           RowBox[{
            SqrtBox["2"], " ", "v", " ", 
            SqrtBox[
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
              "\[Rho]", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", "2"}], " ", "e"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   SuperscriptBox["u", "2"], "+", 
                   SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
               ")"}]}]]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            SuperscriptBox["u", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
           "+", 
           RowBox[{"v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "v"}], " ", 
               RowBox[{"(", 
                RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
              RowBox[{"2", " ", 
               SqrtBox["2"], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ",
                  "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ",", "1"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{"2", " ", "\[Rho]"}], 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
           "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}]]}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{"2", " ", "u", " ", "\[Rho]"}], 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
           "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}]]}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Rho]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}], 
          "+", 
          RowBox[{
           SuperscriptBox["u", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "v", " ", 
              RowBox[{"(", 
               RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
             RowBox[{
              SqrtBox["2"], " ", 
              SqrtBox[
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
                "\[Rho]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "e"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                 ")"}]}]]}]}], ")"}]}], "+", 
          RowBox[{
           SuperscriptBox["v", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "v", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "1"}], "+", 
                SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], 
             "+", 
             RowBox[{
              SqrtBox["2"], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
              SqrtBox[
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
                "\[Rho]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "e"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                 ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"u", "-", "v"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{"u", "+", "v"}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
           "-", 
           RowBox[{
            SqrtBox["2"], " ", "v", " ", 
            SqrtBox[
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
              "\[Rho]", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", "2"}], " ", "e"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   SuperscriptBox["u", "2"], "+", 
                   SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
               ")"}]}]]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            SuperscriptBox["u", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
           "-", 
           RowBox[{"v", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"v", " ", 
               RowBox[{"(", 
                RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
              RowBox[{"2", " ", 
               SqrtBox["2"], " ", 
               SqrtBox[
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ",
                  "\[Rho]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                  ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ",", "1"}], 
    "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513820267393568*^9, 3.513820360153299*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "1", ",", "0", ",", "u"}], "}"}]], "Output",
 CellChangeTimes->{3.513820267393568*^9, 3.513820360356469*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", 
   RowBox[{"v", " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox["u", 
       RowBox[{"2", " ", "v"}]], "-", 
      FractionBox["v", 
       RowBox[{"2", " ", "u"}]]}], ")"}]}], ",", "v", ",", "0"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.513820267393568*^9, 3.513820360690091*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "u", ",", 
   RowBox[{"-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "+", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
            "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}], 
         "+", 
         RowBox[{
          SuperscriptBox["u", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "v", " ", 
             RowBox[{"(", 
              RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
            RowBox[{
             SqrtBox["2"], " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "+", 
         RowBox[{
          SuperscriptBox["v", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", 
               SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], "-", 
            RowBox[{
             SqrtBox["2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"u", "-", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{"u", "+", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "+", 
            RowBox[{
             SqrtBox["2"], " ", "v", " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
          
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["u", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "+", 
            RowBox[{"v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "v"}], " ", 
                RowBox[{"(", 
                 RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
               
               RowBox[{"2", " ", 
                SqrtBox["2"], " ", 
                SqrtBox[
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
                  " ", "\[Rho]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                   ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ")"}]}]}], 
   ",", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{
       RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["u", "2"], "+", 
         SuperscriptBox["v", "2"]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
      RowBox[{
       SqrtBox["2"], " ", "v", " ", 
       SqrtBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
         "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], " ", "e"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
          ")"}]}]]}]}], 
     RowBox[{"2", " ", "\[Rho]"}]]}]}], "}"}]], "Output",
 CellChangeTimes->{3.513820267393568*^9, 3.513820360886827*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "u", ",", 
   RowBox[{"-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "-", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
            "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "4"}], " ", "e", " ", "v", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]"}], 
         "+", 
         RowBox[{
          SuperscriptBox["u", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "v", " ", 
             RowBox[{"(", 
              RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
            RowBox[{
             SqrtBox["2"], " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "+", 
         RowBox[{
          SuperscriptBox["v", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", 
               SuperscriptBox["\[Gamma]", "2"]}], ")"}], " ", "\[Rho]"}], "+", 
            RowBox[{
             SqrtBox["2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "3"}], "+", "\[Gamma]"}], ")"}], " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}]}], ")"}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "e", " ", "\[Gamma]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"u", "-", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{"u", "+", "v"}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "-", 
            RowBox[{
             SqrtBox["2"], " ", "v", " ", 
             SqrtBox[
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
               "\[Rho]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "2"}], " ", "e"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                ")"}]}]]}]}], ")"}]}], "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]", " ", 
          
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["u", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
            "-", 
            RowBox[{"v", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"v", " ", 
                RowBox[{"(", 
                 RowBox[{"1", "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
               
               RowBox[{"2", " ", 
                SqrtBox["2"], " ", 
                SqrtBox[
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", 
                  " ", "\[Rho]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "2"}], " ", "e"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["u", "2"], "+", 
                    SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
                   ")"}]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], ")"}]}]}], 
   ",", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{
       RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["u", "2"], "+", 
         SuperscriptBox["v", "2"]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "-", 
      RowBox[{
       SqrtBox["2"], " ", "v", " ", 
       SqrtBox[
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
         "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], " ", "e"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
          ")"}]}]]}]}], 
     RowBox[{"2", " ", "\[Rho]"}]]}]}], "}"}]], "Output",
 CellChangeTimes->{3.513820267393568*^9, 3.513820361082849*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"ebt", "//", "FullSimplify"}], "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.513820369032012*^9, 3.513820383863572*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "1", "0", "u"},
     {"1", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{"u", "-", "v"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"u", "+", "v"}], ")"}]}], 
       RowBox[{"2", " ", "u"}]], "v", "0"},
     {"1", "u", 
      RowBox[{"v", "-", 
       FractionBox[
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
        RowBox[{
         SqrtBox["2"], " ", "\[Rho]"}]]}], 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "+", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}], 
        RowBox[{"2", " ", "\[Rho]"}]]}]},
     {"1", "u", 
      RowBox[{"v", "+", 
       FractionBox[
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
        RowBox[{
         SqrtBox["2"], " ", "\[Rho]"}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]",
            " ", "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}], 
       RowBox[{"2", " ", "\[Rho]"}]]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.51382040060229*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eb", "=", 
  RowBox[{
   RowBox[{"Transpose", "[", "ebt", "]"}], "//", "FullSimplify"}]}]], "Input",\

 CellChangeTimes->{{3.513820412585874*^9, 3.513820427037327*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "1", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1", ",", 
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{"u", "-", "v"}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"u", "+", "v"}], ")"}]}], 
      RowBox[{"2", " ", "u"}]], ",", "u", ",", "u"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "v", ",", 
     RowBox[{"v", "-", 
      FractionBox[
       SqrtBox[
        RowBox[{
         RowBox[{"-", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
         " ", "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], " ", "e"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
       
       RowBox[{
        SqrtBox["2"], " ", "\[Rho]"}]]}], ",", 
     RowBox[{"v", "+", 
      FractionBox[
       SqrtBox[
        RowBox[{
         RowBox[{"-", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
         " ", "\[Rho]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], " ", "e"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["u", "2"], "+", 
              SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
       
       RowBox[{
        SqrtBox["2"], " ", "\[Rho]"}]]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"u", ",", "0", ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]",
            " ", "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}], 
       RowBox[{"2", " ", "\[Rho]"}]]}], ",", 
     FractionBox[
      RowBox[{
       RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
       
       RowBox[{
        SqrtBox["2"], " ", "v", " ", 
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
           ")"}]}]]}]}], 
      RowBox[{"2", " ", "\[Rho]"}]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.513820428177108*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eb", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.513820451205735*^9, 3.513820455527036*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "1", "1", "1"},
     {"1", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{"u", "-", "v"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"u", "+", "v"}], ")"}]}], 
       RowBox[{"2", " ", "u"}]], "u", "u"},
     {"0", "v", 
      RowBox[{"v", "-", 
       FractionBox[
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
        RowBox[{
         SqrtBox["2"], " ", "\[Rho]"}]]}], 
      RowBox[{"v", "+", 
       FractionBox[
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
        RowBox[{
         SqrtBox["2"], " ", "\[Rho]"}]]}]},
     {"u", "0", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["u", "2"], "+", 
            SuperscriptBox["v", "2"]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], 
         "+", 
         RowBox[{
          SqrtBox["2"], " ", "v", " ", 
          SqrtBox[
           RowBox[{
            RowBox[{"-", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", 
            "\[Gamma]", " ", "\[Rho]", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], " ", "e"}], "+", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 SuperscriptBox["u", "2"], "+", 
                 SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
             ")"}]}]]}]}], 
        RowBox[{"2", " ", "\[Rho]"}]]}], 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "e", " ", "\[Gamma]"}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
        RowBox[{
         SqrtBox["2"], " ", "v", " ", 
         SqrtBox[
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]",
            " ", "\[Rho]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "e"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["u", "2"], "+", 
                SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
            ")"}]}]]}]}], 
       RowBox[{"2", " ", "\[Rho]"}]]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.513820456082701*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"v", "-", 
  FractionBox[
   SqrtBox[
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", " ", 
     "\[Rho]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], " ", "e"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
   RowBox[{
    SqrtBox["2"], " ", "\[Rho]"}]], "+", "a", "-"}]], "Input",
 CellChangeTimes->{{3.513820582618715*^9, 3.513820583276757*^9}, {
  3.51382061967034*^9, 3.513820624955898*^9}}],

Cell[BoxData[
 RowBox[{"v", "+", 
  SqrtBox[
   FractionBox[
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", 
     RowBox[{"(", 
      RowBox[{"e", "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["u", "2"], " ", "\[Rho]"}], "2"], "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["v", "2"], " ", "\[Rho]"}], "2"]}], ")"}]}], 
    "\[Rho]"]], "-", 
  FractionBox[
   SqrtBox[
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", "\[Gamma]"}], ")"}], " ", "\[Gamma]", " ", "\[Rho]", 
     " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], " ", "e"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]], 
   RowBox[{
    SqrtBox["2"], " ", "\[Rho]"}]]}]], "Output",
 CellChangeTimes->{3.513820584032531*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], " ", "e", " ", "\[Gamma]"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["u", "2"], "+", 
          SuperscriptBox["v", "2"]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}], " ", "\[Rho]"}], "+", 
       
       RowBox[{
        SqrtBox["2"], " ", "v", " ", 
        SqrtBox[
         RowBox[{
          RowBox[{"-", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", 
          " ", "\[Rho]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "e"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["u", "2"], "+", 
               SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], 
           ")"}]}]]}]}], 
      RowBox[{"2", " ", "\[Rho]"}]]}], "-", "H"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Rho]", ">", "0"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.513820699204967*^9, 3.513820732874904*^9}, {
  3.513820767785276*^9, 3.513820769304369*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"v", " ", 
    SqrtBox[
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "\[Gamma]"}], ")"}]}], " ", "\[Gamma]", " ",
       "\[Rho]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "e"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["u", "2"], "+", 
           SuperscriptBox["v", "2"]}], ")"}], " ", "\[Rho]"}]}], ")"}]}]]}], 
   RowBox[{
    SqrtBox["2"], " ", "\[Rho]"}]]}]], "Output",
 CellChangeTimes->{
  3.513820704652417*^9, {3.513820744353155*^9, 3.513820770059702*^9}}]
}, Open  ]]
},
WindowSize->{1270, 766},
WindowMargins->{{61, Automatic}, {Automatic, 24}},
FrontEndVersion->"7.0 for Linux x86 (32-bit) (February 25, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 1883, 55, 165, "Input"],
Cell[2453, 79, 405, 6, 31, "Output"],
Cell[2861, 87, 454, 8, 48, "Output"],
Cell[3318, 97, 454, 8, 48, "Output"],
Cell[3775, 107, 405, 6, 31, "Output"],
Cell[4183, 115, 866, 23, 50, "Output"],
Cell[5052, 140, 2158, 62, 60, "Output"],
Cell[7213, 204, 2158, 62, 60, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9408, 271, 1010, 28, 55, "Input"],
Cell[10421, 301, 3811, 117, 91, "Output"],
Cell[14235, 420, 3811, 117, 91, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18083, 542, 3882, 118, 97, "Input"],
Cell[21968, 662, 3743, 117, 91, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25748, 784, 3882, 118, 97, "Input"],
Cell[29633, 904, 3742, 117, 91, "Output"]
}, Open  ]],
Cell[33390, 1024, 96, 1, 32, "Input"],
Cell[CellGroupData[{
Cell[33511, 1029, 1292, 35, 55, "Input"],
Cell[34806, 1066, 2527, 79, 95, "Output"],
Cell[37336, 1147, 2527, 79, 95, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39900, 1231, 2769, 82, 101, "Input"],
Cell[42672, 1315, 2204, 71, 85, "Output"]
}, Open  ]],
Cell[44891, 1389, 571, 18, 32, "Input"],
Cell[CellGroupData[{
Cell[45487, 1411, 2649, 80, 101, "Input"],
Cell[48139, 1493, 2273, 72, 85, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50449, 1570, 908, 26, 77, "Input"],
Cell[51360, 1598, 409, 13, 50, "Output"],
Cell[51772, 1613, 484, 16, 60, "Output"],
Cell[52259, 1631, 491, 16, 68, "Output"]
}, Open  ]],
Cell[52765, 1650, 166, 3, 32, "Input"],
Cell[CellGroupData[{
Cell[52956, 1657, 17157, 494, 584, InheritFromParent],
Cell[70116, 2153, 14604, 421, 333, "Output"],
Cell[84723, 2576, 157, 3, 31, "Output"],
Cell[84883, 2581, 341, 11, 44, "Output"],
Cell[85227, 2594, 7165, 208, 184, "Output"],
Cell[92395, 2804, 7133, 207, 184, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[99565, 3016, 162, 3, 32, "Input"],
Cell[99730, 3021, 4286, 132, 138, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[104053, 3158, 192, 5, 32, "Input"],
Cell[104248, 3165, 3900, 125, 113, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[108185, 3295, 126, 2, 32, "Input"],
Cell[108314, 3299, 4287, 132, 136, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[112638, 3436, 686, 22, 63, "Input"],
Cell[113327, 3460, 977, 34, 70, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[114341, 3499, 1321, 41, 60, "Input"],
Cell[115665, 3542, 684, 23, 62, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
